Time periodicity and dynamical stability in two-boson systems 
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We calculate the period of recurrence of dynamical systems made of two interacting bosons. A 
number of theoretical issues related to this problem are discussed, in particular, the conditions 
for small periodicity. The knowledge gathered in this way is then utilized to propose a notion of 
dynamical stability based on the stability of the period. Dynamical simulations show good agreement 
with the proposed scheme. 



PACS numbers: 

INTRODUCTION 

An important result known as the Recurrence Theo- 
rem P, Q establishes that classical as well as quantum 
systems of constant energy will come back very close to 
their initial state at some (maybe long) time during their 
evolution. A straightforward consequence of the theorem 
is that these systems posses, at least in a broad sense, an 
intrinsic or natural frequency given by the inverse of the 
recurrence's time. Such a frequency characterizes the sys- 
tem's response to external driving. For instance, the am- 
plitude of the response is maximally enhanced when the 
frequency of the driving matches the natural frequency 
of the system This enhancement underlies quite a 
number of important physical phenomena and its under- 
standing is of fundamental interest. Assuming that the 
frequency of the driving is constant, one should minimize 
fluctuations of the natural frequency in order to get an 
amplified response, especially when the parameters of the 
system are subject to small perturbations. The stability 
of quantum dynamics as a result of a small change in the 
parameters has been discussed in several works, for ex- 
ample in chapter 11 of Q or It is seen that in models 
with classical counterparts, the overlap of the quantum 
state of two systems with the same initial condition but 
slightly different parameters remains close to one as long 
as the initial quantum condition is well localized inside 
a stable island of the classical map. The opposite takes 
place when the initial condition is localized in a chaotic 
region. Here we intend to approach the issue of stability 
by studying the behavior of periodicity in a model with 
no direct classical analog. The study contains a moder- 
ate analysis of the period, which lays the foundation for 
the subsequent argument about stability. In absence of a 
classical analog, we base our approach exclusively on the 
cigenencrgies of the Hamiltonian, making little reference 
to the state. We discuss a series of issues in connection 
to the period, particularly, the question of how small the 
period can be as a function of the parameters. The re- 
sults obtained here are relevant in the context of quan- 
tum computation, where some information protocols or 
quantum gates are sensitive to perturbations of the pa- 



rameters and periodic systems can be seen as quantum 
memories Q. 

In the language of the theory of second quantization, 
the quantum state is written with reference to occupation 
modes of unperturbed levels. In this context, let us focus 
on a model featuring two-body interaction such as, 



H = J(oja 2 +q|oi)+y] ^a\a k {a\a k -l)+t k a\a k . (1) 

k=l 

As usual, the exchange term J and the unperturbed 
energies ei and e 2 underline the single body response 
while Ui and U 2 determine the intensity of contact in- 
teraction among particles. The mode operators satisfy 
the usual bosonic relations [&i,a{] = [02, Sj] = 1, etc. 
The unperturbed system can be probed by looking at 
the absorption profile of an incident laser of frequency 
v = I eg — The total number of particles, 



M 



a[ai + a\a 2 , 



(2) 



is a preserved quantity. The proposed system reduces 
to a two-level model when M = 1. Hamiltonian ([1]) can 
be rearranged as, 



H = rja\a\a\di — fia\di + a\a 2 + a 2 ai, 



(3) 



where, for M = 2, the parameters 77, [i and A turn out 
to be, 



3U 2 + C/i A 



2 J 



J 



, A = e 2 -e x . (4) 



In writing the previous identities we have established 
J as our energy unit Q. The state evolution is given by 
(in what follows, we set H = 1), 

|V(*)> =q 1 e- iE i t \E 1 )+q 2 e- iE > t \E 2 )+q 3 e- iE ' t \E 3 ). (5) 
Hence, periodicity arises whenever Q, 



2 



E]T = 27mi, E 2 T = 2nn 2 , E 3 T = 2im 3 . (6) 

Being T the period of the recurrence and ni, n 2 and 
n 3 integer numbers. In this case periodicity is absolute 
as the quantum state recurs identically at regular inter- 
vals and the corresponding evolution operator equals the 
unity operator. Another form of periodicity Q emerges 
by considering instances in which the quantum state re- 
curs up to a phase, i.e., 

m + r))=e-^\m)- (7) 

We refer to this as partial periodicity as the phase fac- 
tor might generate quantum interference effects. As an 
example, let us look at Hamiltonian for a single par- 
ticle, 

A-(ri). («) 

Absolute periodicity occurs whenever the ratio of en- 
ergies is a fractional number, 

E 2 ii+ ^/i_l 2 + 4 n 2 
T can be found from Eq. ([5]), 

T = 2ixn 2 \] — x = 2ir\/—iiin 2 . (10) 

In Eq. © we can define E\ and E 2 so that \n 2 \ ^ |ni| 
and therefore — 1 ^ x < 0. In principle, there is no limit 
on the maximum value of T. Conversely, the minimum 
value is T = 2tt and takes place at // = 0. 

In a similar way, partial periodicity derives from, 



TWO PARTICLES 

Let us now probe these periodicity concepts in a bigger 
system. Hamiltonian Q for M = 2 takes the matrix 
form, 

/ 477-2^ s/2 \ 
H=\ V2 77-M \/2 ■ (13) 
V \/2 / 

The energies are the solutions of the characteristic 
equation, 

E 3 + aE 2 +0E + j = O, (14) 

where, 

a = — 5r] + (15) 

j 8 = 2(2»?-A*)(»7-M)-4, (16) 

7 = 4(277-//). (17) 

From the previous equalities we can see that if an en- 
ergy E is a solution of Eq. ([M)) for a set of parameters 
{?7,/x}, then — is a solution for the set {— i], — /i}- Ad- 
ditionally, given two solutions of Eq. (fl"4"]). £3 and £7, we 
can write, 

£f - £ 3 + a(E% - E 2 ) + f3(E 3 -E) = 0. (18) 
If E3 E, the polynomial can be simplified as, 

E 2 + E(E 3 +a) + E 2 + aE 3 + /3 = 0, (19) 

and the respective solutions provide the unaccounted 
energies, 



E X T 



E 2 t 



2tt, 



(11) Ei = E 3 + a+ y/(E 3 + a) 2 - 4(£| + aE 3 + g 



and therefore, 



2tt 



2- 



Vm 1 



= 2tt 



(1-^ 



(12) 



The first equality is consistent with the view that the 
system's natural frequency is proportional to the differ- 
ence of its two eigenencrgies. Unlike T, r reaches a max- 
imum of t = 7T at /j, = and goes down asymptotically to 
zero as fi — > ±00, where one of the eigenenergies domi- 
nates the spectrum and the Hamiltonian is almost singu- 
lar. This shows that r is maximum when T is minimum 
and that r goes to zero as T goes to infinity. Two level 
systems always display partial periodicity, but not always 
absolute periodicity 



and similarly for E 2 . Absolute periodicity results when 
the energy ratios adhere to, 



Ei ni E 2 n 2 

E 3 n 3 E 3 n 3 



(21) 



One may ask whether given a set of parameters 77 and 
fi the system would display periodicity. This however 
might not be a convenient approach, since in any case we 
can find integers m, n 2 and n 3 showing ratios as close to 
x and y as we want. This is also true for any reasonable 
quantum closed system. Instead, we propose an approach 
in which given a pair of ratios we ask if a set of parameters 
producing these ratios do exist. Following this idea we 
proceed as follows. 
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Standard Procedure 

Let us then consider x and y as known variables. We 
divide Eq. ([20]) by E 3 and get x. y is obtained in a 
similar way. By means of algebraic operations we can 
put the unknown variables in terms of a, 



-a/(l + x + y), 



P 
7 



-pa *, 
■ qa 3 . 



So that p and q are given by, 



(22) 
(23) 
(24) 




p = [x + y + xy)/(l + x + y) 2 , 
q = xyj (1 + x + y) 3 . 

From equations ([To]) and ([IT]) we get, 



7 = g(»7-a)- 



(25) 
(26) 



(27) 



Similarly, using (|2"T1) in combination with ([M]) leads to, 



^1 , 3<1 2 

r? = a I 1 + —a 
Inserting (J25J in (HU) yields, 



(28) 



FIG. 1: Coordinate map indicating the number of set of pa- 
rameters {r/,fj,} delivering energy ratios {x,y}. Red (dark 
gray) region: solutions. Green (light gray) region: 4 so- 
lutions. Blue (continuous) lines: 2 solutions. Yellow dots 
({x = —l,y = 0} and {x — 0,y — —1}): 1 solution. Black 
(dashed) lines highlight the pairs for which there are solu- 
tions with 77 = 0. In these cases, given a couple {x, y} there 
are two solutions with rj = and two solutions with 77 7^ 0. 
The mirror symmetry around y = — x derives from the fact 
that Ei and E2 can be swapped in Eqs. (|21[) . Here we require 



\EA < \E 3 \ and \E 2 



\Es\. Along the edges of the square 



where x or y is — 1, there exists an ambiguity in the choice 
of E3 because two eigenenergies display the same absolute 
value but opposite sign. This implies that the ratio that is 
not —1 can be positive or negative. As a consequence, the 
borders remain identified: {x — — 1, y} ~ {x — —1, —y} and 
{y = -l,x] - {y = -1, -x}. 



= 2a [ 1 + —a 2 
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(29) 



Combining equations ^B, ([23]), j28J and ([29]) we get 
the characteristic equation, 



k 2 t] 6 + 27k 2 i] 4 + ( 243k 2 + ^- ) r/ 2 + 729k 2 = 



q 2 A 3 + 2qA 2 + ApA + 16 = 0, 



(30) 



where A = a 2 . In this way, given a pair of values x 



and y we find p and q from Eqs. (|25[) and ([26]) and in- 
troduce them in Eq. ([3"(J]) . From the solutions we have a 
and therefore the rj and /1 that deliver the energy ratios. 
Strictly, we have 6 solutions, but they come in pairs giv- 
ing energies with opposite sign. In order for a solution to 
be physically acceptable we demand that rj and /i should 
be real. Let us now discuss particular cases where the 
previous method does not apply. 



Particular Cases 



in which we have introduced, 



1 



3 + (x - y) 2 V 3 + {x - y) 2 
Moreover, /j can be found from Eq. (fl~5]) . 



(32) 



Case B. x = 1 or y = 1 (excluding {x = l,y = 0}, 
{x = 0, y = 1} and {x = y = 1}). 



In order to avoid division by zero in Eq. (|18p we modify 
the variables in the following fashion, 



Case A. a = x + y+ l = 0. 



if x = 1 =^ x' = y' 



(33) 



For a pair of values x and y the parameter 77 is a 
solution of the polynomial equation (sec appendix A), 



and analogously when y = 1. In this way the new 
variables admit the standard procedure. 
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FIG. 2: Top: T (arbitrary units) vs. x, keeping y constant. 
Bottom: T vs. y, keeping x constant. In both cases the grid 
slice is 1/97. As can be seen, the smallest values of T are 
found close to the origins. 



Case C. {x = 0, y = 1} or {x = 1, y = 0}. 



= r\. Non- vanishing 
and x (or y) = 1 



Then q = 7 = 0, j3 — —4 and a 
energy values satisfy E 2 + r\E — 4 = 
is a solution only if 77 = ±4i. 

Case D. {x = y = 1}. 



The spectrum of the Hamiltonian is three fold degen- 
erate. Hence the eigenvalue equation is of the form 
(E 3 — E) 3 = 0. Comparing with Eq. (fT4|) we infer 
-ZE = a,3E 2 = /3 and -E 3 = 7. It then follows 
E 6 + 6E 4 + 12E 2 + 16 = 0, which has no real solutions 
and therefore no real- valued 77 and /1 produce a fully de- 
generate spectrum. 

Fig. Q] shows a map classifying the coordinate space 
according to the number of valid solutions encountered 
for every pair of ratios {x, y}. The maximum number 
of solutions is 4, usually coming from 2 real solutions of 
Eq. (f30)) . Along the negative side of the axes the Hamil- 
tonian becomes reducible with /1 = 2rj, hence only two 
solutions are possible. The two instances with 1 solution 
correspond to 77 = /j, = 0. 



Recurrences 

The condition for absolute periodicity is given in Eq. 
([5]) while partial periodicity can be determined through, 



E 2 t = 



2irN 2 , E 3 t = c/) + 2itN 3 , (34) 



in such a way that N 2 and N 3 are both integers. If 
T = iVr, being an integer, we can infer from Eqs. ((BJ) 
and ([3"i]). 




FIG. 3: Partial period (arbitrary units) against x for M = 2 
and M = 1 (Eq. (fT2|) - dashed blue line). We have inten- 
tionally restrained the vertical axis to < r < 10. Interact- 
ing Hamiltonians can display longer r than non-interacting 
Hamiltonians. For this particular graph we used grid slices 
ranging from 1 to 1/500 in order to scan the xy square of Fig 
[T](see appendix B). The inclusion of finer grids would produce 
a denser, but not very different, graph. This is because finer 
grids are probabilistically related to bigger periods. 



27mi 
N 



, N 2 



n 2 - m 
N ; 



n 3 - n x 
N 



(35) 



In order to find T and r, we first determine 774, n 2 
and 113 from the rationals x and y in such a way that 
there is no common divisor greater than 1 among the 
three generating integers. Simultaneously, the values x 
and y are used to encounter corresponding cigencncrgics 
following the previously discussed methodology. We are 
then entitled to employ any of the identities of Eqs. © 
(in particular we employ E 3 T = 2irn 3 ) to get T using 
one of the eigenenergies. Finally, the integer N results 
as the greatest common divisor of n 2 — n\ and n 3 — n\. 
As can be seen, large values of the integer n 3 are more 
likely to produce large T and r. 

Fig. [5] presents sample diagrams of T against x and 
y. A complete depiction would be quite more intricate. 
We point out that T decreases as x and y approach zero. 
This behavior appears to be generic, laying minimum 
values of T around (or in) the origins. Following T across 
the line {x = —l/n,y = 1/ti} for n = 5,6,7, etc. we 
find it goes down asymptotically toward T — ny/2 as 
11 — > 00. The same minimum value can be analytically 
found at {x = — 1/2, y — 0} and its equivalent. Both 
instances suggest a relation between minimum T and the 
Hamiltonian matrix being or becoming singular. 

Fig. [3] pictures the intricate relation between r and x 
obtained by testing energy ratios over the xy square of 
Fig. [T]as explained in appendix B. Cooperative systems 
may display longer partial periods than non-interacting 
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FIG. 4: Blue lines: set of parameters for which the period 
is stable under small changes in the direction indicated by 
the red lines. Inset, r (arbitrary units) against the smallest 
of \x\ and \y\. It follows that as r goes to zero one of the 
eigenenergies becomes small in comparison to the others. 



systems, therefore periodicity can help identify interact- 
ing phases. This behavior is most likely due to the capac- 
ity of many-body systems to develop complex dynamics 
on account of the increased number of effective states. 
Also of interest is the chaotic aspect of Fig. [3j meaning 
that small intervals in xy space do not necessarily map 
onto small intervals in r space. This lack of continuity, 
which also characterizes T, takes place because neigh- 
boring values of x and y do not always correspond to 
neighboring values of rii, ni and 713. For instance, {x = 
1/3, y = 2/3} is not far from {x = 97/300, y = 201/300}, 
but the sets of integers that generate each pair are widely 
different. Recall that the period depends directly on the 
integers. In addition, there can be common divisors be- 
tween 7J2 — ni and 713 — n\ and this might further affect 
the continuity of r. Finally, it can be seen from the inset 
of Fig. |4]that when r becomes small at least one of the 
ratios x or y approaches zero, once again suggesting that 
in a Hamiltonian displaying small periodicity one of the 
eigenenergies is much smaller than the others. 

STABILITY 

As introduced here, periodicity, either total or partial, 
is a characteristic of the Hamiltonian. It does not de- 
pend (except when the initial state is an eigenstate of 
the Hamiltonian) on the state. This encourage us to 
ask whether one can change the Hamiltonian parameters 
without affecting the period. This is what happens for 
one particle when fi = 0, as can be seen from Eq. (jT2"j) : 



For two particles, we can start by arguing that the pe- 
riod is stable whenever both x and y arc stable under 
changes of the parameters. Moreover, we can consider 
the energy ratios as functions of the parameters, x(ij,fj,) 
and y(?7,/i), in such a way that maximum variations oc- 
cur in the direction of the function gradients and zero 
variations take place in the directions perpendicular to 
the gradients. In this way, for any set of parameters rj 
and /i one can always identify a direction of change of 
the parameters along which one of the ratios is stable. It 
then follows that in order for T to be stable the direc- 
tions of zero change of both ratios must coincide, i.e., the 
gradients must be parallel: 

Vi + XVy = 0, (37) 

where Vx = f^?7+ an d similarly for Vy. if and 
p, are unitary vectors in the parameter space and A is 
a real number. This problem is equivalent to finding 
the extreme values of x(rj, fi) subject to the condition 
2/(77, /z)=constant, or the other way around. In this con- 
text A takes the roll of a Lagrange multiplier and the anal- 
ogy applies as long as the involved functions are smooth. 
The extremes of x can be worked out from Fig. [TJ It is 
found, however, that not all extreme values imply paral- 
lel gradients. For instance, along the edges, where x or 
y is —1, the identification of borders on each side of the 
axes generates a discontinuity in the first derivative of 
the ratios. In addition, for the extremes located around 
the neighborhood of the axes at least one of the param- 
eters diverge toward infinity. Let us remark that there 
is no solution of Eq. (|30l) along the positive side of the 
axes. In general, we found the gradients are parallel only 
along the blue lines of Fig. Q] located in between green 
and red regions. 

Fig. U shows curves indicating the parameters as well 
as the directions of change corresponding to zero varia- 
tion of the period. Far from the origin the curves ap- 
proach (but do not seem to touch) straight lines given 
by simple expressions. In one case the direction of zero 
change aligns with the direction of the curve as we get 
away from the origin. In the other case the direction of 
zero change becomes constant with an angle of inclina- 
tion 9 satisfying tan w 1/ \/2. 

As can be seen from Fig. [5j when the parameters are 
shifted in the direction of zero change the dynamics of 
corresponding systems look similar. Conversely, when 
the change of parameters occurs in the direction of the 
gradient, the dynamics soon diverge. This behavior is 
consistent with our study and shows that in some cases 
the period alone characterizes the global profile of the 
system evolution. 

We have studied time periodicity as well stability prop- 
erties of interacting bosonic systems. We stressed the 
difference between the usual notion of periodicity, where 
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Time 



FIG. 5: Mean number of particles (a\ai) against time (ar- 
bitrary units) for three slightly different cases. Top. r\ — 
1.878615 and fx = 2.939702. Middle. n = 1.778615, fi = 
2.775950. Bottom. 77 = 1.878615 and /i = 2.714882. The 
middle-graph corresponds to a set of parameters in the curve 
of stability of Fig. 2] Top- and bottom-graphs show the evo- 
lution of the system when such parameters are shifted in the 
direction of zero change and the direction of the gradient, 
respectively. 



the ratios of the Hamiltonian energies become commen- 
surate, and partial periodicity, where the state recurs up 
to a phase factor. Both forms of periodicity were ex- 
plicitly established for one and two particles in a model 
described by two bosonic modes. The results suggest 
a connection between minimum periodicity and the fact 
that one of the eigenvalues becomes small against the 
other eigenvalues, especially for the two- particle case [Toj | . 
Likewise, we observed that the stability of the period docs 
not only depend on the parameters but also on the direc- 
tion of change of such parameters. Hence, we presented 
a diagram showing the set of parameters as well as the 
directions of change for which the period stays constant. 
These results were consistent with simulations of the dy- 
namics. 

An issue of the method previously applied to find r 
is the assumption T = Nt for the two-particle case. It 
may be possible to relax this assumption and derive r 
using only Eq. (f3~4")l . This would in principle lead to a 
richer stability diagram including parameters for which 
t is stable but T is not, e.g., r\ = fx = 0. Similarly, we 
feel that the method presented here is robust, although 
slightly dependent on the dimension of the Hamiltonian. 



Appendix A: Derivation of Eq. (1311) . 

When a = the solutions of Eq. (|19p can be written 

as, 



-E 3 - y/-3£f - 4/3 -E 3 + y/-3E* - 4/3 

2 ' E2 ~ 2 

(38) 

Then it follows, 



Replacing ([39]) in ([13]) with a — gives, 

k/3 3/2 = 7, (40) 

where k is given by Eq. (|32[) . Likewise, \i = |?/\ By 
replacing such a value of \i in Eqs. (fTtlj) and (|T7|) we find, 

P = -4 (j + l) , 7 = \n- (41) 

Substituting these values in Eq. (j40|) and after some 
algebraic steps we arrive at Eq. (|3"Tj) . 

Appendix B: Generation of commensurate ratios. 

In order to generate pairs of ratios compatible with 
the coordinate square of Fig. [TJ we first establish an 
integer 77,3 = 1,2,3,..., which is also the inverse of the 
grid slice. For a given 713, we check —713 < m < — 1 
and n\ < n% < n$ — 1. The values of 713 are introduced 
into a computer routine in increasing order starting from 
713 = 1. Subsequently, n 1 and are introduced. The 
integers are used to get x and y which are in turn used 
to check for Hamiltonians with corresponding eigenen- 
ergies. As the coordinate square is symmetric, we scan 
only the section x < 0. If for a given triplet {rti, 77,2, 713} 
there exists a maximum common divisor greater than 1, 
the triplet is discarded. Notice that in this way the scan 
of the coordinate square is performed using several su- 
perimposed grids, avoiding repetition of pairs {x,y}. 
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